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The nanoparticle orientation in fluid systems can be correlated with the rotational diffusion and is widely 
used to tune the physical properties of functional materials. In the current work, the controllability of the 
orientation of a single rigid carbon nanotube in a fluid is investigated by imposing a linear shear flow. 
Molecular dynamics simulations reveal three forms of anomalous behavior: (i) "Aligned orientation" when 
the nanotube oscillates around a particular direction which is close to the flow direction at a small angle of 
about 10° in the velocity-gradient plane; (ii) "Interrupted orientation" when the oscillation is interrupted by 
a 360° rotation now and then; (iii) "Random orientation" when 360° rotations dominate with the rotational 
direction coinciding with the local fluid flow direction. The orientation order is a function of the Peclet 
number (Pe). The results show that the correlation between Pe and the orientation order from the 
two-dimensional model does not apply to the three-dimensional cases, perhaps due to some anomalous 
behavior and cross-section effects. This work provides clear pictures of the nanoparticle movement that can 
be used to guide particle manipulation techniques. 

Scientists are devoting more and more efforts to understanding the dynamics of active nanocomponents in 
soft matter systems for nonequilibrium conditions. Single-particle tracking experiments^"^ have enabled 
real-time monitoring and have provided detailed information about the particle diffusion and directed 
motion. The rotary Brownian-type motion of a non-spherical nanoparticle, such as a carbon nanotube (CNT), 
which has been referred to as rotational diffusion in contrast to translational diffusion, has drawn considerable 
attention recently^"^°. An external force field provides a confining potential which restricts the random nano- 
particle rotation^ The current study analyzes the combined effect of an applied moment and the Brownian 
torque that drives a single CNT to fall into some intriguing patterns, one of which is the orientation type. 
Orientation control of an elongated particle can benefit the in vivo particle transport^^ or biosensor designs^^. 
The present study focuses on nanoparticles in large systems with this analysis of single particle orientation 
transferable to alignments of particles in dilute solutions where inter-particle interactions are not significant. 
CNT based nanocomposites and nanofluids have many useful features, with enhanced mechanicaP^'^^, elec- 
tricaP^'^^, opticaP^ and thermal properties'^. However, in nanofluids'^''', the nanotubes arrange themselves 
randomly, which limits further improvements of the thermal conductivity of such nanofluids. The CNT's 
anisotropy suggests that an aligned distribution may further enhance the physical properties which will facilitate 
the processing and designing of advanced nanocomposites and related materials'^ '^. 

For carbon nanotube suspensions, the general constraint for the dilute regime where a study can focus on a 
single particle is cL^ < 1, where c is the number density and L is the mean length of the nanotubes'^. Furthermore, 
treating the CNT as a rigid rod-shaped object allows direct comparisons with classical theories. Jeffrey in 1922'^ 
gave a theoretical treatment of a rigid ellipsoid in a sheared viscous fluid which showed that, in the absence of 
Brownian motion, the ellipsoid would rotate indefinitely in a single parameter family of closed orbits, known as 
Jeffery orbits. A series of studies have discussed the Jeffery orbits'^"'^ and the nontrivial effects of the very weak 
rotary Brownian motion on the rheological properties of suspensions^^ or single particle movements^\ Molecular 
dynamics (MD) studies^' have further confirmed Jeffery's theory by studying the motion of a low aspect ratio 
CNT in simple shear argon. This study showed that Brownian motion does account for the slightly different 
rotational time periods compared with the theory. Tannous^^ used a different method to derive the two dimen- 
sional (on the velocity-gradient plane) probability density functions from Langevin simulations for rigid solid 
objects. Tannous observed a preferred oriented distribution of rigid solid objects around the flow direction and 
the orientation order could be expressed as function of the Peclet number {Pe = y/Dr, where y is the shear rate and 
Dr is the rotational diffusion coefficient). This conclusion was confirmed for liquid flows inside a planar channel 
geometry^^ 
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Figure 1 | (a) Schematic diagrams of the initial configuration of the simulation system and (b) definition of the orthogonal coordinates, x, y, z, 
longitudinal angle, ^, and latitudinal angle, in the simple linear shear flow; (c) The angle range of the longitudinal angle, ^, in the x-j plane. The flow is 
along the x direction, while y is the gradient direction. 



The studies in the literature stifl do not clarify under what circum- 
stances a rigid sohd object or a CNT would exhibit periodic rotations 
or would locate along an oriented direction^^'^^. This lack of under- 
standing is the result of two problems. One is the lack of a clear 
physical view of the three-dimensional (3D) movement of the object, 
which cannot be obtained from two-dimensional theoretical treat- 
ments, whfle the other is due to insufficient observations and ana- 
lyses of how different physical attributes affect the orientation. 
Actually, there is stfll uncertainty about whether the definition of 

stifl applies to 3D rotary motion of an object. This work addresses 
these issues by investigating the orientation of a single rigid nanotube 
in a fluid undergoing a simple linear shear flow with the focus on the 
3D orientation and the quantitative influence of various factors. This 
study uses a low aspect ratio CNT which can be treated as a rigid body 
to allow direct comparisons with theoretical results. The results with 
this simplified model can also be compared to experimental results 
with very large aspect ratio CNTs to analyze the effect of the rota- 
tional diffusion coefficient^'^. Molecular dynamics simulations are 
used for their simplicity in modeling the atom- atom interactions, 
the direct monitoring of the nanotube's Brownian motion and the 
abflity to provide vivid pictures of how the nanotube actually moves. 
Unlike existing works, this study investigates the nontrivial coupling 
between the rotation, the orientation and the random Brownian 
motion of a nanoparticle, which results in various types of motion 
not observed in earlier studies. 

Results 

The MD simulation system was set up using orthogonal Cartesian 
coordinates, x, y, z, with longitudinal and latitudinal angles cp and 0 



defined as in Figures 1(a), (b). The superimposed simple linear shear 
fiow is iflustrated in Figure 1(b) and given by Eq. (1), 

Vx=yy.Vy = Vz = ^. (1) 

where y stands for the shear rate. The carbon nanotube is surrounded 
by fluid argon and its long axis is initially parallel to the y direction. 
(In Supporting Information, the initial configuration is changed and 
no influences on the results can be observed: Figure SI and Table SI) 
Although the simulations were performed in a 3D frame, the discus- 
sion mostly focuses on the x-y plane, which means that the dir- 
ectional angle is actually the longitudinal angle (p. Considering the 
fore-aft symmetry of the capped nanotube, the angle range is defined 
as shown in Figure 1(c). When the nanotube is paraflel to the flow 
direction x, (p is equal to 0°. In the velocity gradient direction, (p may 
have a discontinuous jump from +90° to —90° and is negative for 
nanotube "A" and positive for nanotube "B" as shown in Figure 1(c). 
The base case is defined here with various factors then changed 
individually to examine their influences. The capped nanotube had 
the armchair (5, 5) configuration with a base case diameter, d, of 
0.7 nm and length, L, of 6.8 nm. The density, p, of the argon was 
1763 kg/m^ and the temperature, T, was 300 K, which corresponds 
to a supercritical state. The simulation domain had a size of L^X LyX 
Lz = 10.6 X 10.6 X 10.6 nml The shear rate was equal to 2 X 
10^° s"^ for the nonequilibrium simulations. 

The orientation type of the nanotube will be affected by the prop- 
erties of the nanotube, the fluid and the shear rate. To simplify the 
discussion, the length, L, of the nanotube is changed first to study the 
motion with the other factors held constant at the values in the base 
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Figure 2 | Time-varying instantaneous angle ^ for (a) 1 = 8.1 nm; (b) 1 = 1.8 nm; (c) L = 0.9 nm. Here, L denotes the length of the nanotube. 



SCIENTIFIC REPORTS | 4:6120 | DOI: 1 0.1 038/srep061 20 



2 




Figure 3 | Time-varying instantaneous angle ^ with a range of 0° to 360° plotted in polar coordinates for (a) 1 = 8.1 nm; (b) 1 = 1.8 nm; (c) L = 

0.9 nm. Here, L denotes the length of the nanotube. Here, the angle range is 0° to 360° to display a full circle of rotation, different from the definition in 
Figure 1(c). 0° means that the nanotube is parallel to the y direction. 



case. Three types of motion appeared for the time-varying instant- 
aneous angle, (p, as shown in Figures 2(a) -(c). Firstly, for L equal to 
8.1 nm, the nanotube oscillates around a particular angle, approxi- 
mately + 10° from the flow direction (Figure 2(a)). Secondly, for L 
equal to 1.8 nm, most of the time, the nanotube still has a preferred 
orientation angle of around +10° with several circular rotations 
mixed into the oscillations now and then (Figure 2(b)). Finally, the 
nanotube with the smallest length rotated continuously all the time 
with no fixed angle (Figure 2(c)). As nanotube becomes shorter, the 
orientation of the nanotube becomes less ordered and the nanotube's 
orientation type can be classified into three forms as "aligned ori- 
entation", "interrupted orientation" and "random orientation". This 
anomalous behavior is interesting since there are no existing theories 
to predict this behavior, especially for the unusual interrupted ori- 
entation. The aligned orientation state is of particular interest due to 
its relevance to particle alignment. These observations can be better 
understood by plotting the instantaneous angles for the three forms 
in polar coordinates as in Figures 3(a)-(c). Now, the angle range is 0° 
to 360° to display a full circle of rotation. (Here, 0° means that the 
nanotube is parallel to the y direction which differs from the angle 
definition in Fig. 1(c) which applies to the rest of the manuscript.) 
The 360° rotations of the nanotube are always clockwise as seen in 
Figures 3(b) -(c), as was observed earlier^^'^^. This rotational direction 
can be anticipated for the superimposed linear shear given by Eq. (1) 
which creates a local clockwise vorticity while the viscous drag by the 
surrounding fluid induces a torque on the nanotube. In addition, no 
clear boundary exists between the interrupted and random orienta- 
tions. The transition from the interrupted state to the random ori- 
entation state occurs when the nanotube's rotations dominate over 
the oscillations. Actually, these observations suggest that the occur- 
rences of the orientation types are not always consistent and the 
complex combinations of mixed orientation and rotation types can- 
not be predicted theoretically^l 



The underlying mechanisms for these three types of orientation 
were further investigated by examining the nanotube's three dimen- 
sional movements as shown in Figures 4 and 5. Figures 4 (a) -(c) show 
the projection of the time- varying ratio of the nanotube's length on 
the x-y plane while Figures 5 (a) -(c) show the time- varying angular 
velocity, co^. At all times, the projection of the ratio is close to unity 
for the longest nanotube whose orientation type is aligned orienta- 
tion (Figure 4(a)). This means that the nanotube is almost in the x-y 
plane with a very small latitude angle 0. For the random orientation 
case with L = 0.9 nm, the projection of the ratio varies greatly 
between 0 and 1 (Figure 4(c)). Thus, a nanotube that rotates out of 
the x-y plane with a larger angle 0, will be more likely to experience a 
360° rotation away from the preferred orientation angle. As the 
latitudinal angle, 0, becomes larger, the nanotube rotates faster 
around the z axis, co^ is approximately ±0.06 rad/ps for L = 
8.1 nm (Figure 5(a)) and ±0.6 rad/ps for L = 0.9 nm 
(Figure 5(c)). Moreover, for the random orientation type, the time 
period for angle (p is approximately 834 ps. For an object rotating 
around Jeffery orbits with negligible Brownian effect^^, the time per- 

2n f 1 \ 

iod is about 4490 ps as estimated from T = — [ Vg-l — ) where y is 



the shear rate and = 1 .24 - 



y 



is the equivalent aspect ratio^^ 



Thus, these results indicate faster rotations for our simulations, 
which is reasonable considering the small friction imposed on the 
nanotube from the surrounding argon in the supercritical state. The 
interrupted orientation type is like a combination of the other two 
types. Specifically, seven rotations can be differentiated from the 
dominate oscillations (Figure 2(b)) with a marked decrease in pro- 
jection of the ratio (Figure 4(b)) and a sudden distinct increase in the 
absolute value of co^ as shown in Figure 5(b) with only one exception 
(No.®). As the projection ratio is equal to cosO, the latitudinal angle 
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Figure 4 | Time-varying projection of the ratio of the nanotube's length on the x-7 plane for (a) 1 = 8.1 nm; (b) 1 = 1.8 nm; (c) L = 0.9 nm. Here, L 
denotes the length of the nanotube. 
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Figure 5 | Time-varying angular velocity, ^y^, for (a) 1 = 8.1 nm; (b) 1 = 1.8 nm; (c) L = 0.9 nm. Here, L denotes the length of the nanotube. 



(0° ~ 90°) is generally less than 40° (the ratio is larger than 0.8) for 
the oscillations and well above 40° when the rotations occur for these 
three cases. Thus, to have a good orientation, the nanotube should be 
confined to the x-y plane with a relatively small latitudinal angle. 
Otherwise, the 360° rotation is more likely to develop with a larger Wz 
than during the oscillations. 

The orientation variations upon changing the diameter, d, of the 
nanotube or the density, p, or temperature, T, of the argon are quite 
similar to that of changing L, with the results presented in the 
Supporting Information (Figures S2-S7). The qualitative trend is 
that increasing d or T, or decreasing p relative to the base case will 
cause the nanotube orientation type to shift from aligned to inter- 
rupted and then to the random orientation type. The corresponding 
3D variations of the rotational direction, projection ratio and angular 
velocity are all quite similar. In addition, the nonequilibrium system 
approaches equilibrium as the shear rate decreases. At high shear, 
occasional oscillations around a particular angle can still be distin- 
guished and the 360° rotation is always clockwise for the random 
orientation type (Figure 3(c)). However, for very low shear, the rota- 
tion is random and unpredictable (Figure 6). This phenomenon 
arises due to the increasing relative strength of the random 
Brownian couples relative to the mechanical torque due to the shear 
flow near equilibrium. 



The Peclet number. 



Pe = y/Dr, 



(2) 



in which y is the shear rate and is the rotational diffusion coef- 
ficient, indicates the relative strength of the hydrodynamic interac- 
tions and Brownian forces and is, thus, a crucial parameter indicating 
the orientation order. To exhibit the results in terms of Pe, the rota- 
tional diffusion coefficient, Dj., must first be defined. Cao and Dong^^ 
calculated this diffusivity using equilibrium and nonequilibrium 
molecular dynamics methods. A detailed discussion of how to cal- 
culate Dj. is presented in the Supporting Information (Figure S8 and 
Table S2). Pe can then be obtained from Eq. 2 for the simulation 
cases. The angle probability density distributions, P{(p), has isolated 
peaks for the two different lengths in Figure 7 indicating that the 
nanotube shows preferred orientations around these particular 
angles. Higher peaks and narrower distribution ranges result in a 
more ordered orientation type for larger Pe. The higher peak appears 
for the larger Pe and the probabilities for angles away from this peak 
are nearly zero. For smaller Pe, the probabilities are almost the same 
over the entire angle range, indicating no preferred orientation angle. 
The present MD results do not agree well with the results of 
Tannous^^ shown in Figure 7. Although PMD(^max) = 0.629 and 
P[33](^niax) " 0.608 are quite close for Pe = 4.36, the large difference 
between PMD(^max) = 2.93 and P[33](^max) = 1-33 for Pe = 49.3 
should be noticed. 
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Figure 6 | Time- varying instantaneous angle ^ with a range of 0° to 360° 
plotted in polar coordinates for shear rate y = 1 X 10^ s~\ 
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Figure 7 | Angle probability density distributions for two different 
nanotube lengths of 0.9 nm and 3.1 nm compared with theoretical 
results^^. 
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Figure 8 | (a) Dependence of the orientation order parameters T^x and /I on the aspect ratio. The cases with different lengths, L, or diameter, dy are 
identified by different symbols; (b) Dependence of the orientation angle, ^, on the aspect ratio. The inset is a sketch of the nanotube placement to show 
that the nanotube has a positive angle with respect to the flow direction, as represented by the solid outline rather than the dotted outline; (c) Dependence 
of the orientation order parameters, T^xy on the Peclet number compared with theoretical results^^. The cases of changes in the other factors with a fixed 
diameter are differentiated from the cases with different diameter with the other factors fixed by different symbols. 



The nanotube orientation type can be further defined by the order 

tensor^^'^^: 



(3) 



where N is the total number of nanotubes and is the unit axial 
direction vector. The subscripts / and; are selected between x andy 
since the orientation is only considered in the x-y plane. Thus, T is a 2 
X 2 symmetric matrix. If the diagonal element Ta = 0.25, the time- 
varying angles exhibit a completely random distribution. If Ta = 1, 
the nanotube is perfectly oriented along the / axis. The largest eigen- 
value of T, symbohzed as /I, is actually the orientation order along the 
most preferred angle, i.e. orientation angle (p^, which is the corres- 
ponding eigenvector. The orientation order parameters, T^x, and the 
largest eigenvalue X of T were calculated for various aspect ratios 
(Figure 8(a)). The armchair configuration of (5, 5) to (12, 12) cor- 
responding to a diameter change of 0.7 nm to 1.7 nm was used with 
the length varied from 0.9 nm to 8.1 nm. In all cases, 1 is larger than 
Txx, indicating that the preferred direction is not parallel to the x axis. 
Both of these parameters become closer to unity as the aspect ratio 
increases. In Figure 8(b), the orientation angles show no clear 
dependence on the aspect ratio, with all these angles in the interval 
of (+2°, +18°), approximately consistent with Tannous^^ Positive 
orientation angle mean that the nanotube is more likely in the state 
shown in the solid sketch in the inset of Figure 8(b) rather than the 
state shown by the dotted outline. Detailed results for the orientation 
order parameter and orientation angle are shown in the Supporting 
Information (Figure S9) for various densities, temperatures and fluid 
shear rates. The results show that with increasing fluid density, 
decreasing temperature, or increasing shear rate, the nanotube tends 
to become more oriented along a specified direction. The results also 
show that the orientation angle has a very slight dependence on these 
factors and is about +10° relative to the flow direction. When 1 is 
larger than about 0.95, the orientation type has the aligned orienta- 
tion; when 0.6<>^<0.95, the orientation type has the interrupted 
orientation; and when 1 is less than 0.6, the orientation type is the 
random orientation. 

The orientation order parameter, T^x, was then calculated for 
various Pe for a large number of cases with the results for different 
nanotube diameters indicated by different symbols in Figure 8(c). 
One of the main conclusions in Tannous^^ was that a positive cor- 
relation exists between Pe and the orientation order parameter, T^x- 
As can be seen, the MD results in Figure 8(c) clearly differ from 
Tannous's work with no positive correlation between T^x and Pe. 
In particular, the diameter changes produce opposing trends. In 



Figure 8(c), d is varied from 0.7 nm to 1.7 nm for two nanotube 
lengths of 4.3 nm and 6.8 nm. The diameter increase leads to a larger 
Dj. as shown in Table SI (Supporting Information) and a smaller 
orientation order as shown in Figure 8(a). This phenomenon may 
challenge whether the expression for Pe or its relationship with the 
orientation order established in Tannous^^ are valid. 

Thus, in the same shear flow field, a more easily rotating nanotube 
with a larger Dy. will not necessarily lead to a less ordered time- 
varying orientation with significant differences between the orienta- 
tion states and the theoretical treatment in Tannous^^. In fact, the 2D 
study by Tannous showed that the rigid solid object motion and the 
object's cross-sectional area are both negligible in the plane of the 
flow. As the nanotube is free to rotate in the 3D frame in the present 
study, so its movements may become quite complicated, as shown by 
the three anomalous orientation types discussed here. The additional 
effect due to the cross section may also be responsible for the differ- 
ences from the theory, particularly as evidenced by the differences in 
the results as the diameter changes. 

Discussion 

To conclude, this study reveals anomalous orientation variations of a 
rigid carbon nanotube in a sheared fluid in molecular dynamics 
simulations. First, the orientation variations can be classified into 
three forms as the "aligned orientation", "interrupted orientation" 
and "random orientation", corresponding to the decreasing orienta- 
tion order in each type. The relatively small projection of the nano- 
tube length ratio onto the x-y plane means that the 360° rotation can 
more easily occur with a relatively larger angular velocity, co^, which 
will interrupt the oriented state even leading to the state of no ori- 
entation at all. Second, the nanotube has a preferred time-varying 
directional angle, i.e. a longitudinal angle, ^, of around +10° from 
the flow direction. The time-varying orientation type is found to be 
affected by the length and diameter of the nanotube, as well as the 
density, temperature and shear rate of the fluid. Third, the Peclet 
number serves as an important scaling parameter for processing the 
results. The results show that the Peclet number does not positively 
correlate with the orientation order, especially as reflected by the 
results showing that the diameter increase does not reduce the order 
parameter. This breakdown may arise from the non-trivial cross 
sectional effect of the diameter and the nanotube's complex 3D rota- 
tional movements during the three types of motion. The nanoparticle 
rotational patterns could potentially be exploited to control the 
nanoparticle movement in bio-transport processes, the fabrication 
of materials with exceptional properties and other applications. We 
believe this work will benefit the study of dilute CNT suspensions 
and can be extended to other nanoparticle suspensions focusing 
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mainly on analysis of single particle motion. The CNT is assumed to 
be rigid as in many experimental investigations to allow direct com- 
parisons with existing theories of rotational diffusion. Although the 
current results cannot be directly applied to concentrated suspen- 
sions, these results will serve as a useful point of reference for theor- 
etical formalisms and experiments with concentrated suspensions 
that must incorporate the correct physics that holds in this limit. 
The anomalous behavior observed here encourages us to devote 
more efforts into modifying and developing theories to describe these 
behavior types, which will in turn serve as the foundation for more 
up-to-date design and analysis techniques. 

Methods 

The rodlike CNT nanotube molecule is capped at both ends to prevent argon atoms 
from entering the nanotube. The MD package LAMMPS was used for the calcula- 
tions^*^. The forces between the argon atoms and between the argon and carbon atoms 
were calculated based on the Lennard-Jones (LJ) pair potential: 



The basic parameters were tiAr-Ar = 0.341 nm, SAr-Ar = 1-67 X 10"^^ J, tJc-Ar = 

0. 357 nm, and Sc-Ar ^ 1-97 X 10"^^ The C-C interactions were not considered 
since the nanotube was assumed to be a rigid body^^, so a larger time step of 10 fs 
could be used in the calculations. The time-consuming calculations of the inter- 
particle interactions were reduced by using a cutoff distance of 0.77 nm. The NVT 
canonical ensemble was used with the Nose-Hoover thermostat. 

The superimposed simple linear shear flow was realized by utilizing the SLLOD 
algorithm^\ The Lee-Edwards boundary condition"^^ was applied in the y direction, 
which is the velocity gradient direction with periodic boundary conditions for the 
other two directions. The excess heat produced by the shear flow was extracted by the 
thermostat. The thermostat was only applied in two directions {y, z) and not in the 
flow direction {xY^'"^^. The results were collected after 3 X 10^ time steps (3000 ps) to 
allow time for the flow to develop. 
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